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Abstract 

We use the constraints arising from primordial nucleosynthesis to bound a putative electric 
charge density \e\ng of the universe. We find \nq/n.y\ ^ 10"^'^, four orders of magnitude more 
stringent than previous limits. We also work out the bounds on Uq in models with a photon mass, 
that allows to have a charge density without large-scale electric fields. 
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One of the fundamental parameters of the universe is its electric charge. It is usually 
assumed to be exactly vanishing, but it is of course desirable to have observational evidences 
of the hypothesis of a neutral universe. In this paper we will discuss constraints on Ug, defined 
in such a way that \e\nq is the electric charge density of the universe (e is the electron charge). 
It is convenient to normalize this quantity to the photon number density, so we define the 
parameter 

V,^^ (1) 

with 

% = ^ (2) 

where T is the temperature, and ^(3) ^ 1.20. 

Speculations about a possible electric charge of the universe go back to the work of 
Lyttleton and Bondi |^ , who showed that an excess charge could account for the expansion 
of the universe. They assumed that the excess charge was due to a tiny charge imbalance 
among the proton and the electron, qp + qe ^ 0. Laboratory limits 0, Q on this imbalance 
rule out the assumption of Lyttleton and Bondi. Their argument can be somehow inverted 
demanding that the gravitational attraction among cosmological objects is larger than 
the electromagnetic repulsion leads to the bound 

1^,1 ^ 10-^^ (3) 

More recent speculations involving an excess charge are exotic suggestions of violation of 
electric charge conservation. Grand unified theories imply violations of baryon and lepton 
numbers, but conservation of electric charge is usually believed to be exact as a result of 
gauge invariance. Nonetheless one has the possibility of a (necessary small) breakdown of 
electric charge conservation at high energies Also, we may have an excess charge in 
theories with extra dimensions, first introduced by Kaluza and Klein and subject to a 
recent strong revival 0. The extra dimensions should be compactified to a small size, since 
our world appears four dimensional. To probe them one can consider the universe in its 
first instants; the energy is so high that one has a truly higher-than-four dimensional space. 
Subsequently, the universe evolves and cools down but some electric charge could remain in 
the ordinary space. Some authors have built classes of models where there is charge leakage 
into infinite extra dimensions [HI. 
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It is therefore interesting to examine the imphcations of a charge density riq and to 
find constraints on it. Orito and Yoshimura |^ examined some of these imphcations. They 
considered the anisotropies generated by the large-scale electric field that would be produced 
by a net charge density. Their best bound is obtained using observational limits on cosmic 
ray anisotropies: 

1^,1 ^ 10-'' (4) 

In this letter we show that a much more stringent bound can be obtained by considering 
the primordial nucleosynthesis period of the early universe and the change in the primordial 
helium production that would be originated by a charge density present at the nucleosyn- 
thesis era. 

The physics involved in the primordial nucleosynthesis period is well understood and the 



theoretical predictions of the primordial yields of light elements are robust 0. The 
agreement with observation is considered one of the pillars of modern cosmology. For our 
purposes, it is useful to recall here the simple arguments that allow to understand the main 
physical features of helium production |P, |ll|. In the early universe, at time t <^ 1 s and 
temperature T ^ 1 MeV, neutrons and protons are in kinetic and chemical equilibrium 
due to weak interactions. At a lower T ^ Tf, the rate of these interactions becomes less 
than the expansion rate of the universe and they go out of equilibrium. Then, the relative 
neutron-to-proton density is given approximately by 



n 

.P. 



-Ql^f ^ 0.24 (5) 



/ 



where 



Q = mn-nip= 1.29 MeV (6) 

We have taken Tf 0.9 MeV. After freeze-out, there is a decrease in the neutron number 
due to /5-decay, 

n ^ p e~ u (7) 
until the time tne where all the helium is produced. Then 

Here we have introduced the neutron lifetime Pl t„ = 887 s and we have taken tffe ~ 400 s. 
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Since nearly all neutrons are processed into helium, the expected mass fraction Y of ^He 

is 

2 (n/p) Hp , ^ 

Y ~ — ^ 0.26 (9) 

1 + in/p)He ' 



The actual theoretical prediction for Y has to be obtained by a numerical code [12, 13] that 
solves the relevant set of ordinary differential equations, and indeed gives a result not far 
from d^). 

Primordial nucleosynthesis is known to be a useful tool to constrain non-standard physics 
01 . We now apply it to the issue of a charged universe. A charged particle in such universe 
has a potential V (with the condition ^ = when nq = 0) and thus it has a momentum p 
related to its total energy E by 

{E - Vy =p^ + m^ (10) 
With the introduction of an effective mass m*, the dispersion relation can be written as 

E^=p^ + {m*f (11) 
We now assume V <^ E and that the charged particle is non-relativistic. We get 

m* ~ m + (12) 



A similar phenomenon is at the basis of the MSW effect |14]; eqs(|T0|-12) in the relativistic 



case are discussed in [|T5 



The horizon distance an is defined as the propagation time of light since t = 0, and then 
for r > an the charge density could not interact with the charged particle. Thus, to the 
potential V contributes all the charge inside a sphere of radius an 

V = ± r (drAnr^)^ (13) 
JO r 

To summarize, a net charge density \e\nq induces effective masses for the charged particles 
present in the universe at that time 

nip m* = nip + 5m 

me —>■ m*^ = me — 6m (14) 
(positrons are affected with the opposite sign than electrons), where from ([1^) and (|13D 

5m = 271 aa\n^rjq (15) 



As we are in a radiation dominated universe we will put an = 'it. 

A non-vanishing 6m alters the predictions of primordial nucleosynthesis for the helium 
yield, and this leads to a bound on 6m and hence on rjg. Before presenting our numerical 
results let us show that we can get an approximate expression for the change in Y due to 
the 6m shift, using the same type of simple analysis that we used in our discussion of the 
standard Y, from eq. (j^) to (|^). We work at first order in 6m. 

There are two main sources of change in Y. The first is through the dependence of Y in 
eq. (H) on Q at freeze-out, 

5Y 1 6m 6m 

~Y ~ 1 + {nlv)f tJ ^ MeV 



(16) 



where we have introduced the numerical values for Ty and {n/p)f. 

The second is due to the dependence of r„, the lifetime of process (^, on the proton and 
electron masses. When decaying into a proton of mass m* and an electron of mass m*, the 
neutron lifetime (for zero temperature) depends on masses as 



— oc m* X{q 



(17) 



/q* 
dxx {x — q*Y Vx^ — 1 



60 ^/'^ 



l(2g**-9g*^-8) 



q* \n{Jq*^ -1 + q*) 



(18) 



with q* = {mn — m*)/ml 

After some algebra one finds the change in r„ 



5Tn 



5-(g-l) 



'1 dX 
X dq* 



6m 
m. 



(19) 



with q 



{mn — mp)/me. From (H) it is easy to see that a shift 5r„ induces a change 6Y 

6Y 1 



- — 5Tn 

{n/p)He 



(20) 



Y 1 

Introducing the numerical values for the parameters that appear in (|19]) and (|20|) , we end 
up with 

Am 

(21) 



6Y 
~Y 



MeV 



To find the expressions (16) and (21) we have made the assumption that 6m is indepen- 
dent of time, which is not true, since an = 2t and ~ T'^ ~ t^^/"^ in ([T5|). Still, we are able 



to find the right order of magnitude for our bound on rjg if we proceed as follows. Since ([I6D 
is approximately valid at freeze-out, we evaluate 6m in (|T5|) using t ^ 1 s, T ^ 1 MeV and 
get 

6m ^ 10^^ r]g MeV (22) 

Similarly, (^) is approximately valid during the period of neutron decay. Now we evaluate 
6m in ([TsD with typical values for this period, t ~ 100 s, T ^ 0.1 MeV, and get 

6m ^ 10^2 r]g MeV (23) 

Comparing (ESI) and (PSf), we see that it is the change in Y due to (5r„, eq. (P5|), that 



will dominate the whole effect. The origin of this dominance is that, as we see from (P!3[), 

increases with time. To estimate the order of magnitude of 
the bound on r^g, we put (|23D in (^) and allow, for instance, a 10% change in Y . We get 
\r,g\ < 10-43 

To find a precise bound, we have implemented the modification (0) and ([T5| ) in Kawano's 



version 



jl^ of the primordial nucleosynthesis code of Wagoner to obtain the effect of 
a density Ug on the helium abundance. The abundance depends on the laboratory neutron 
lifetime r„, the number of neutrinos A^^, and the baryon-to-photon ratio rjB- We use r„ = 
887 ±2 s 0, A^",^ = 3 0, and for rjB the value extracted from CMB anisotropy measurements 
p!6| . Unfortunately, this value is still not completely settled but the hope is that in the 
future it will be with more refined experiments. For the time being, we take 

10"^° <riB< 10"^ (24) 



a generous range that embraces the observations [|16 . 

The prediction for F as a function of rjg is shown in Fig.|I], for the two extremes of the 
range (|2^). (The experimental error in r„ introduces a negligible change in our prediction.) 
The linear dependence is due to the expected dominance of the first-order in rjg. Our bound 
is stringent because Y is measured with relatively small error. The observational status is 



discussed in We adopt the range 

0.228 <Y < 0.248 (25) 



that the authors of [10| claim is "95% C.L.". As we see from Fig.||, combining the experi- 



mental range (p5|) with the predictions for Y in the range (^4]) constrains the charge density 



of the universe. The bound is 

- 1.1 X 10"^^ <Vq< 0-8 X 10"^^ (26) 

We notice that it is four orders of magnitude better than and that it is not far from our 
estimate using (^31) . 

A non- vanishing rjg would also affect the yields of the other light elements D, '^He, and 
^Li. However, since these yields are not as well measured as ^He, taking them into account 
could not significantly improve our bound (pB]). 

Even if an hypothetical charge density of the universe Ug has to be so tiny, one may worry 
about the induced large-scale electric fields that appear when 7^ 0, no matter how small. 
In fact, there are models having 7^ but no large-scale fields. This can be achieved 
by endowing the photon with a small mass and thus having an electromagnetic 



interaction of finite range A. Bounds, like (^), that are based on the effects of a large-scale 
electric field are very much weakened in this kind of models (to get @) the authors of @] 
consider a scale ~ 500 pc for the electric field). 



Let us show that our bound ([26|) , when 7^ 0, is only modified by one order of magni- 
tude. First, we notice that the mass shift 6m ~ t^/^ shown in (|15]) is valid when an = 2t < \. 
However, for later times, 6m gets contribution only up to a radius A, 

<5m = — e(3)r/,AT3~r3/2 (27) 

71 

It follows that 6m increases until t = X/2 and afterwards it decreases. 

Since helium production finishes when t ^ tne, it is clear that our bound (^Gf ) is still valid 
for A > tne- For smaller A, we expect to find bounds on that are less severe. Obviously, 
for A — the bound would disappear. However, A is subject to the experimental constraint 

A = — >10^m (28) 



coming from studies on torques on a toroid balance |T8[. Similar limits are obtained from 
measurements of the Jovian magnetic field [|l^. We notice that the lower limit A = 10^ m is 
about the same as)f:/~ls = 3x 10^ m. Thus, to estimate the bound on \rig\ we can use ( ^2]) 
since it is valid for t ^ 1 s. We expect to get \rig\ < 10"'^^, so that in models with a photon 
mass we expect to have a bound that is only about one order of magnitude worse than (p6D. 
We have modified Kawano's code with the mass shift 6m shown in (27) to find the 



precise change in Y. By demanding that the predicted Y is in the experimental interval (p5D, 



and allowing for t]b itl the range (0), we are able to limit rjq. In Fig.^, we show our bound 
as a function of A and also the experimental limit (pBI). We confirm that for A ^ 10^^ m we 
get our previous limit ( pB]) and we also see that the constraints relax for smaller A. We can 
put a limit on 1]^ in models with a photon mass, which clearly is the one for A = 10^ m. It 
reads 

- 1.6 X 10"^^ <Vq< 1-8 X 10"^^ (29) 
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0.26 ^ 




FIG. 1: Predicted value of y as a function of 7]q = Ug/riy, for r]B = 10~^ (left-upper line) and 
for rjB = 10^^'^ (right-lower line). The two horizontal lines are the observational limits (p5|). The 
allowed range of Y is not shadowed. 




FIG. 2: Upper and lower bounds on rjg when electromagnetic interactions have a range A = 

The experimental limit on A (^) is displayed as a vertical line. Non-excluded values are in the 

non-shadowed part of the figure 
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